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Abstract 

We present a simple way of separating the overlap between the soft and collinear 
factorization formulae of QCD squared matrix elements. We check its validity ex- 
plicitly for single and double unresolved emissions of tree-level processes. The new 
method makes possible the definition of helicity-dependent subtraction terms for 
regularizing the real contributions in computing radiative corrections to QCD jet 
cross sections. This implies application of Monte Carlo helicity summation in com- 
puting higher order corrections. 
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1 Introduction 



The majority of available techniques for the computation of radiative corrections in Quan- 
tum Chromodynamics (QCD) relies on the universal factorization formulae of QCD matrix 
elements for the emission of unresolved (soft and/or collinear) partons, which describe the 

^ On leave from University of Debrecen and Institute of Nuclear Research of the Hungarian 
Academy of Sciences, Hungary. 
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universal structure of the infrared singularities [1,2]. These factorization formulae have 
played an essential role in devising completely general algorithms [3-8] for the cancellation 
of infrared singularities when combining the tree-level and one-loop contributions in the 
evaluation of jet cross sections at the next-to-leading order (NLO) accuracy. More recently 
a new algorithm has been proposed [9] that can also be extended to the computation of 
the next order (NNLO) radiative corrections [10,11]. The essence of these algorithms is the 
definition of suitable approximate matrix elements that have the same infrared singular 
behaviour as the matrix elements themselves. The factorization formulae however, cannot 
be directly used as such approximate matrix elements for two reasons. On the one hand 
the factorization formulae are unambiguously defined in the strict unresolved limits. On 
the other hand the soft and collinear factorization formulae overlap in regions of the phase 
space where the unresolved parton is simultaneously soft and collinear to a hard parton. 
The naive application of the factorization formulae as subtraction terms leads to double 
subtractions and as a result to uncancelled infrared singularities. 

There are two ways to cope with the double subtractions. One can either introduce compen- 
sating terms for the double subtractions as for example in the classic paper [12], or define 
single subtractions that smoothly interpolate between the soft and collinear regions, as in 
the dipole scheme [8]. Both techniques have also been extended to NNLO computations. In 
Refs. [13-15] antennae subtraction terms have been introduced that smoothly interpolate 
among all doubly unresolved regions, while the systematic way of accounting for the overlap 
among the various factorization formulae for double unresolved emissions has been worked 
out in Ref. [16]. Thus the problem of double subtractions could in principle be considered 
solved. Yet in this letter we present a third solution to this problem, which is very simple 
and can be extended to any order in perturbation theory easily. It also offers some algorith- 
mic advantages if one aims at automatizing the computation of QCD radiative corrections. 
Furthermore, it allows the introduction of helicity dependent subtraction terms as we show 
in this paper. 



2 Separation of the soft and collinear singularities at NLO 

The separation of the soft and collinear singularities can be obtained from the different 
physical picture of the two cases. In a physical gauge the collinear singularities are due 
to the collinear splitting of an external parton [17]. The overall colour structure of the 
event does not change, the splitting is entirely described by the Altarelli-Parisi functions 
which are a product of colour factor^ and a kinematical function describing the collinear 
kinematics of the splitting. The emission of a soft gluon is just the opposite. It does not 

^ Eigenvalues of the quadratic Casimir operators of the emitted collinear gluon in the represen- 
tation of the parent parton, or the normalization factor of the colour charges in the case of gluon 
splitting into fermions. 
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Fig. 1. General structure of infrared factorization at any perturbative order. (From Ref. [18] with 
permission of S. Catani.) 

affect the kinematics (momenta and spins) of the radiating partons, but it affects their 
colour because it always carries away some colour charge. As a result it leads to colour 
correlations that can be percieved as a soft gluon cloud around the event. 

The analytic expression for the soft factorization of the QCD matrix elements can be 
derived using the soft-gluon insertion technique [2,18]. To describe this technique, we use 
the colour-state notation [8] for the amplitudes and the operator notation [16] for taking 
the soft limit of the amplitude. In the colour-state notation the amplitude for n external 
partons is represented by the ket vector \M.^{pi, . . . ,Pn))- In the limit, when momentum 
p!f becomes soft, the m + 1 parton amplitude fulfills the following factorization formula (for 
the precise meaning of the operator Sr and also that of the operator Cir used below, we 
refer to Ref. [16]) 

Sr{cr\M^^,iPr, ...))= e^{pr)J,{r, e)\Mj. . . )) , (2.1) 

where Cr is the colour-index of parton r. Both the amplitude and the soft current Jr can 
be expanded in perturbation theory. 



\M) = IM'^'^) + \M^'^) + . . . , 



Ja{r,e) = g^fi' 



(2.2) 
(2.3) 



In this letter we concentrate on the leading terms. The leading-order contribution to the 
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soft current is 

4°nr) = En^, (2.4) 

where = (c^lT^'' is the colour-charge operator of parton k if the emitted gluon has 
colour index Cr- The amplitude is a colour- singlet state, therefore colour conservation can 
be expressed as 

m 

J2Tl\M^ip„...,p^)) = 0, (2.5) 
fc=i 

which implies conserved eikonal current, 

m 

P'rJf{r)\M^{p,, . . . = E Tl\M^{p,, . . .,pj) = . (2.6) 

k=l 



The soft limit of the squared matrix element can be obtained by squaring Eq. fl2.ll) . 

mm 1 

SAM^:,UPr, ■■■)\' = Sna^fi'^ E E E eMKiPr)^Si':ir)\M'^Z,^,^{. ..)\\ (2.7) 

4 = 1 fc=l pol. 

with 

Stk\r) = ^, (2.8) 

where Sir = '2'Pi-Pr and \-M.m- (i k)\'^ denotes the colour-connected squared matrix element. 
In evaluating the square, we use axial gauge with a light-like gauge vector n'^ {n^ = 0) to 
sum over the gluon polarizations, which leads to the gluon polarization tensor d'^'^: 



d>^''{Pr,n)=J2e,{Pr)e:{Pr) = -9'" + ^^-—^. (2.9) 



pol. Pr-^ 



As discussed above, in a physical gauge the collinear singularities are due to the collinear 
splitting of an external parton, which remains true in the soft limit. The physical picture 
of collinear and soft emissions suggests that the collinear part of the soft gluon emission 
can be singled out by the diagonal term in the double sum {i = k) in the soft limit 
of the squared amplitude. As mentioned, colour-conservation implies conserved eikonal 
current (see Eq. (12.61) ). therefore, the gauge terms in the gluon polarization tensor do not 
contribute. However, if we want to separate the diagonal terms from the rest, we must keep 
the gauge dependent terms. A short calculation yields our new formula for the soft limit of 
the squared matrix element: 
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S.|A<S^i(p.,...)r = -8vras/x'^E 

i=l 



Sir Sr 



(2.10) 



where Sin = 2pi-n and S,-i.{r) = g^^S'i'^{r), i.e. 



S.k{r) = ^. (2.11) 

Exploiting colour conservation, the gauge- dependent terms cancel and we obtain the well- 
known form of soft-gluon factorization. It is easy to check that the terms in the first line, 
containing the colour correlations, are finite when parton r is coUinear to either i or k, 
while those in the second line become 

C,.S.|AlSii(p„ . . . )|^ = 87ra.f,'^Y.T'. -\M^^\. ■ .)\\ (2.12) 

in agreement with Ref. [16]. Therefore, the soft-coUinear contribution has been separated 
in the second line. 

In order to define the soft limit in Eq. fl2.10p explicitly, we have to fix the gauge vector n^. 
The choice for the gauge vector defines the infrared subtraction scheme employed for the 
computation of radiative corrections because if we want to avoid the double counting of the 
soft-collinear contribution, we have to be able to identify exactly the same expression for 
the soft-collinear contribution both in the soft limit (colour- diagonal piece in Eq. (12.101) ) 
and in the collinear limit expressed in terms of momentum fractions (like in Eq. (I2.12p ). 
The momentum fractions in the Sudakov parametrization of the collinear limit are defined 
only in the strict collinear limit, and in order to define subtraction terms one has to extend 
their definition over the whole phase space. This can be done using a reference momentum 
as 

_ SiP _ SrP , - 

Zi — , Zy — , yz.LO) 

S{ir)P S{ir)P 

where S(jr)p = SiP + Srp so that Zi + Zr = I- If one aims at setting up a subtraction scheme 
valid at any order in perturbation theory, has to be chosen such that S(jr)p must not 
vanish in any of the multiple unresolved (soft and/or collinear) regions of the phase space. 
In Ref. [9] was chosen to be the total incoming momentum Q'^ of the event, when 
Zi/zr = Siq/srQ, SO Eq. f l2.12p reads 

"I 9 e 

C,rSr\M^:UPr. ■■■)\' = STra^/i^^E^?- — I-^^^H- • • )r • (2.14) 

i—l Sir SrQ 

This expression would be identical to the colour-diagonal piece in Eq. fl2.10p if the gauge 
vector were chosen ra^ = Q^, which is not possible in light-cone gauge. We may however. 
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choose n'^ as 

n'' = ariQ^ - brp^) , (2.15) 

such that n"^ = and Sm = SrQ = 2pr-Q. These two requirements determine = 1 and 
K = <5^/ SrQ uniquely. This choice is equivalent to the Coulomb gauge in the center of mass 
frame (rest frame of Q^). With this gauge vector and using colour-conservation, we can 
rewrite Eq. fl2.10p as 



S.|-MSi-i(Pr,...)r = -8vr«s/i'^E 



i=l 



2si 



ik\ 



2sfeQ l|A^(o) 



! 2 Si 



■\M 



(0), 
m > 



(2.16) 



that is, formally we can make the simple replacement Q^. Therefore, in order to avoid 

double subtractions, we simply drop the colour-diagonal terms from the soft factorization 
formula. 



Sr\M 



(0) 

m+1 \Fr^ ■ ■ ■ 



—Sirasfi 



2e 



m -I m 
i=l k^i 



2siQ 



2SkQ 
^rQ^kr 



\M 



(0) , 

m;{i,k) ^ 



(2.17) 

Using this purely soft limit, the collinear limit and the phase-space factorizations of Ref. 
[9], we arrive at the same NLO subtraction scheme as in Ref. [9]. Thus it appears that 
our proposal for separating the soft and collinear limits has not brought any advantage 
as compared to the usual technique of subtracting both and adding back a proper soft- 
collinear compensation term. Note however, that the separation of the soft and collinear 
limits based on the colour structure makes the procedure very simple to any order in 
perturbation theory. Furthermore, it allows for defining helicity-dependent subtractions 
and consequently, Monte Carlo summation of helicities in the computation of radiative 
corrections to jet cross sections. We discuss these two points in turn. 



3 Separation of soft and collinear singularities in multiple infrared emissions 



The overlapping structure of soft and collinear divergences in multiple infrared emissions 
becomes very complex rapidly with increasing number of unresolved partons. Already for 
two unresolved partons there are triply overlapping regions. The separation of these requires 
very careful analysis of the various infrared limits [16] and is a rather laborous excercise. The 
simple rules defined in the previous section can be applied to automate this cumbersome 
procedure as we discuss below. 

We consider first the doubly soft-coUinear limit when momenta and are collinear and 
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the gluon momentum p^, (s 7^ i, r) is soft0 The factorization formula may be written in 
the form [18] 



2e\2 



Sir 



(3.1: 



CSir-s\M^^\2{Pi,Pr,Ps, • • = 

X (A^(^Hp(..), • • OI^ [J^IM ^)d,.iPs, n)J1„){s, e)] Pj^X^,, z,, k^, e)\M^^\p^,r), ■ ■ •)) 



where the soft current J^.^^(s, e) is given by Eqs. (12.31) and (12 ■4p with the replacement r ^ s. 
The subscript {ir) serves simply to remind us that in the summation in Eq. (12.41) . k may 
take the value {ir), in which case the summand is 



j^s P(ir)^i _ ^j^s _|_ rps^ iPi Pr ) fj, 



(ir) 



P{ir) -Ps 



{Pi+Pr)-Ps 



(3.2) 



Retracing the steps leading to Eq. (12.161) . in particular using colour conservation to make 
the replacement n^^ Q'^, we find 



,Pr,Ps 



2t\2 



2s 



kQ 



^ir j—if^-/.j \ ^sQ^js ^sQ^ks 



2 s 



x(A^^Hp(,,), . . .)\PfX{z.. zr, kr, e)T,T,\M^:^^{p^„), • • •)) 
(X(^)(P(,,), . . .)\Pf;i{z,, zr, k^; e)\M^^\pi„), . . .)) 



rrtl 
- J- {ir) 



2 s 



{ir)Q 



{M^^\p{.r), . . ■)\P}:iM, zr, k^- e)\MZ\p^ 



,(0) 



(0)/ 



•)) 



.(3.3) 



It is straightfroward that the product in the first two lines is finite when p^ is coUinear to 
any other momentum that appears in the matrix elements on the right hand side, including 
PI Furthermore, defining the momentum fractions as in Ref. [10], we see that the third 
and fourth lines just reproduce the double and triple collinear limits of the doubly soft- 
collinear factorization formula (Eqs. (4.36) and (4.32) in Ref. [16]) respectively. 

Next, we discuss the case of double soft parton emission. The soft-gluon insertion rules are 
applicable in any order of perturbation theory [18], therefore, the soft-factorization formula 
for the amplitude can easily be given. For instance, for two soft partons the soft current 



We note that if parton s is a fermion, then the squared matrix element does not have a leading 
(doubly-unresolved) singularity in this limit. 

^ We remind the reader that unless explicitly indicated, we include parton (ir) in the summations. 
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has the expansion 



where the leading-order contribution in hght-cone gauge is [18] 



(3.4) 



•r Pipi. rps Pji- 
j 



Pi-Pr ' Pj-P. 

m / 



r rjis Piy PifJ. _|_ rps rpr Pif^ Pi' 
' i i 



Pi-Ps Pi-P{rs) 
SI Pia 



Pi-Pr Pi-Pirs), 



i Pi-P{rs) Pr-Ps 

™ T) 1 

~ Pi-P{rs) Pr-Ps 



In Eqs. (133]) and ^ pf^^^ = p^f + and 

1 



V|3f^u{Pr,Ps 



{Pr - Ps)f39^LU + Psfi9uf3 - Pry g lip 



(3.5) 

(3.6) 



(3.7) 



In order to separate the soft and collinear terms, we apply the same procedure as in the 
case of single unresolved emission, which means taking the square of the soft current and 
separating the colour-diagonal terms. In the case of soft-gg emission the algebra is relatively 
simple. Using (Ti-Tr)^ {Tk-Tr) = T^Ti-Tk, we obtain 



i k 



Pi-Pirs) 



Pk-P(rs) 



where n"|(p,,,Ps) is the quark contribution to the discontinuity of the gluon propagator. 



, .^yPrP^'s+PsPr-r'^Pr-Ps]. 



(3.9) 



Separating the colour-diagonal contributions we find 



(0), 



y [ y r d^aiPirs),n) ^^p^^ ^ df3^{p^rs),n)pl 



(3.10) 
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It is now straighforward to check that in the hmit when the momenta p^, and either 
or are colhnear, the first hne in Eq. fl3.1UI) does not contain leading colhnear singularities, 
while the expression in the second line becomes 



which agrees with the expression in Eq. (4.37) of Ref. [16]. 
We choose the gauge vector similarly as in Eq. fl2.15p . 

Requiring = 0, Sm = SrQ and Ssn = Ssq (so that also S(rs) 
1 



R 



= M (1 _ 



c. = ^(l-i?) 

Sr.e: 



-\- Zg Si{j's'^S'psyZf ~h Zg 



(3.11) 



(3.12) 

n = ■^(rs)^), we can determine 



R 



\ 



1 



SrQSsQ 



. (3.13) 



We can now use colour-conservation to see that the formal substitution can again 

be applied to obtain 



m ( m 



'-'-qq \Pr^ Ps) 
Pi-P(rs) Pk-P(rs) 



'. 2 



.(3.14) 



The soft-collinear term is now separated in the term proportional to T^. This term can 
simply be dropped because it also appears in the Altarelli-Parisi splitting function of the 
q qiq'rO's colhnear splitting [18] if the momentum fractions are defined as in Ref. [10]. 

The same procedure works also when two soft gluons are emitted. In this case it is more 
convenient to rewrite the double-current in an equivalent form in terms of colour-charge 
anticommutators and commutators: 



= EE om,T,^}^;f.(pr,P.) + E[n,T|]c;,(p.,p. 



(3.15) 



where 
and 



A%{Pr.Ps 



Pifj, Pju 



Pi-Pr Pj-Ps 

Piad°''^{p(rs),n)V(3f,u{Pr,Ps) 1 (Pr " Ps) Pif^Piu 



Pr-Ps Pi-Pirs) 



^Pi-iPr + Ps) Pi-PrPi-Ps 



(3.16) 



(3.17) 
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In evaluating the square, we use the following colour identities: 



Ti-T^Sii — Ti-Ti {6ik + Ski) 



5jk 



+ U,{T^ T\ Tl 6u - Tt Tl 5,,) 



[Tl,TtnT^,,T^]=CAT,.T, 



(3.18) 

(3.19) 
(3.20) 



The terms proportional to the structure constants fate in Eq- (13.181) are antisymmetric, 
while the kinematic factors that these terms multiply are symmetric when i with k and j 
with / are simultaneously interchanged. Therefore, these do not contribute to the square. 
The remaining terms contain four-fold and two-fold summations. That with four sum- 
mations has single diagonal {J^iT'lY.jHi^jT j-Ti) and double diagonal {Y.iTfY.j^i'T'^ or 
Y.iTfT'^) terms, which separate the soft-collinear (Eq. (4.45) in Ref. [16]), doubly-collinear 
(Eq. (4.43) in Ref. [16]) and abelian triply-collinear (Eq. (4.40) in Ref. [16]) contributions 
in the doubly-soft emissions, respectively. The contribution with two summations has single 
diagonal terms (I]jT^), which separate the non-abelian triply-collinear pieces (Eq. (4.41) 
in Ref. [16]). Again, the colour non-diagonal pieces do not have leading singularities in the 
various collinear limits. We choose the gauge vector as in Eq. fl3.12p . 



4 Monte Carlo treatment of helicity summation in NLO computations 



In Refs. [19,20] a Monte Carlo integration over a phase variable for reproducing the helicity 
sums in the squared matrix element was introduced in order to save CPU time when 
computing multijet cross sections. While this approach has been found useful for computing 
cross sections at the LO accuracy, its extension to NLO computations is hampered by the 
explicit summation over the helicities in the subtraction terms used in any of the known 
NLO calculations. With the separation of singularities presented here, the Monte Carlo 
treatment of the helicity summation becomes possible by keeping the helicity states for the 
unresolved partons. 

The Monte Carlo treatment of the helicity summation requires helicity-dependent subtrac- 
tion terms. Although the soft and collinear limits of helicity amplitudes are well-known [21], 
the overlap between these cannot be determined at the amplitude level. In taking the square 
of the amplitude, keeping the helicity-dependence in the collinear subtractions is straight- 
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forward. In order to keep the helicity-dependence in the soft subtractions, one needs soft 
terms from which the colhnear singularities are subtracted in a hehcity-independent way. 
Our recipe does precisely that: the soft-coUinear terms are identified in the colour sum. In 
order to define hehcity-dependent soft subtractions, we do not substitute the polarization 
tensor d^jj for the summations over gluon polarizations. The collinear contributions can 
nevertheless be separated in the colour-diagonal terms as discussed in Sect.|2l 

We start from the soft factorization formula Eq. (12.71) . but now we do not sum over the 
polarizatons of gluon r or the helicities of the other m partons. Separating the colour- 
diagonal terms as before, we have 



S.|-M2^i(p.\ 



2=1 



m -1 



(4.1) 



Next, notice that the double summation over i and k in the first line of Eq. ( 14. ip effec- 
tively symmetrizes S^l^{r) in its Lorentz indices, while Sil'{r) in the second line is already 
symmetric in /x and u. Thus only the symmetric part of e^{j)r^n)e^^{jpr,n) contributes. 
However 

^)^v)(Pr, ^) = \ {<^^iPr, n)e^^{pr, Ti) + e^{pr, n)e^^{pr, n)) = ^df,^{pr, n) , (4.2) 
and so by exactly the same steps that lead to Eq. (I2.16p . we obtain 



Sr|A4m+i(p^, . . . 



i=l 



m 



2SiQ 2SkQ \ I ^ (0) ^ 



SrQSir ^rQ^kr 



\M 



m;{i,k)\' ' ' 



Sir SrQ 



■\M 



(0)/ 



(4.3) 



Eq. (14.31) shows that the soft factorization formula is independent of the polarization of 
the emitted soft gluon. Performing the summation over the helecity A of the soft gluon 
(a multiplication by two) as well as the helicities of the rest of the partons, we trivially 
recover the expression in Eg. f|2.16l) . Dropping the soft-coUinear term in the second line, 
we find the helicity-dependenil I purely soft subtraction term, that clearly does not contain 
leading singularities when parton r is collinear to either parton i or parton k. The helicity- 
dependent collinear subtraction terms, including the soft-coUinear contributions, can be 



^ Recall that in Eq. (j4.3p . the helicities of all partons (not just the helicity A of the soft gluon) 
are fixed and not summed over. 
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obtained by squaring the collinear factorization formulae for helicity-amplitudes [21,22]. 
Helicity-dependent subtraction terms may also be defined using the antennae factorization 
expressions of Refs. [23,24]. 

We note that we have checked the validity of Eq. (14. 3 p explicitly for the processes e~^e~ 
qqg and e~^e~ QQgg starting from the expressions for the relevant helicity amplitudes as 
given in Ref. [25]. 



5 Summary 

We have defined a new method for separating the soft and the collinear singularities in 
the QCD factorization formulae. The rules of the method are very simple. One starts 
with the soft-gluon insertion rules for finding the soft limit of the amplitudes. In taking 
the square of those we do not exploit colour-conservation for cancelling the gauge terms 
that appear in the physical polarisations of the soft gluon, rather we separate the colour- 
diagonal contributions. The colour non-diagonal contributions are free of leading collinear 
singularities, while the collinear limit of the diagonal contributions lead to the known 
singular expressions. 

For the gauge vector we may choose a light-like vector whose space-like component 
points into opposite to that of the unresolved gluon r and require that Sm = s^q, which 
ensures that we can perform the formal substitution Q^, where is the total 

four-momentum of the event. This amounts to using Coulomb gauge. We have shown that 
this gauge can be easily generalized to any order in perturbation theory. Choosing the 
momentum fractions in the collinear subtractions as in Ref. [9], the terms separated in the 
colour-diagonal contributions can be identified also in the collinear subtractions, therefore, 
can be dropped. Other choices for the gauge vectors are also possible, but we do not discuss 
that further in this letter. 

This technique can be automatized easily and applied in any order of perturbation theory. 
It also facilitates the use of Monte Carlo helicity summation in the computation of the 
radiative corrections, therefore, can lead to significant reduction of CPU time when there 
are many partons in the final state, which is the most interesting case for the LHC. We 
should mention that a Monte Carlo treatment of colour summation also results in significant 
reduction of CPU time. Such treatment of colour in NLO computations is also facilitated 
by this new method. 

Our method could be useful also for improving parton shower Monte Carlo algorithms. The 
algorithms that are implemented presently treat the colour correlations approximately in 
the large N^. limit and sum up the leading and the next-to-leading logarithms at leading- 
colour accuracy. The summation of the subleading logarithms with exact colour treatment 
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is rather difficult. However, one can go beyond the leading-colour approximation systemat- 
ically by considering the subleading colour contributions pertubatively instead of exponen- 
tiating them. To do this we have to introduce subtraction terms for the subleading colour 
contributions in the parton shower algorithm. The technique presented here can be used to 
define splitting kernels and the corresponding counterterms for a parton shower algorithm 
since the separation of the singularities is governed by the colour structure, which provides 
a good control over the large logarithms and the colour structure simultaneously. We shall 
elaborate these ideas in separate publications. 
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